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ABSTRACT 

We present new wide-field photometry and spectroscopy of the globular clusters (GCs) around 
NGC 4649 (M60), the third brightest galaxy in the Virgo cluster. Imaging of NGC 4649 
was assembled from a recently-obtained HST/ACS mosaic, and new Subaru/Suprime-Cam 
and archival CFHT/MegaCam data. About 1200 sources were followed up spectroscopically 
using combined observations from three multi-object spectrographs: Keck/DEIMOS, Gem¬ 
ini/GMOS and MMT/Hectospec. We confirm 431 unique GCs belonging to NGC 4649, a 
factor of 3.5 larger than previous datasets and with a factor of 3 improvement in velocity 
precision. We confirm significant GC colour bimodality and find that the red GCs are more 
centrally concentrated, while the blue GCs are more spatially extended. We infer negative 
GC colour gradients in the innermost 20 kpc and flat gradients out to large radii. Rotation 
is detected along the galaxy major axis for all tracers: blue GCs, red GCs, galaxy stars and 
planetary nebulae. We compare the observed properties of NGC 4649 with galaxy formation 
models. We find that formation via a major merger between two gas-poor galaxies, followed 
by satellite accretion, can consistently reproduce the observations of NGC 4649 at different 
radii. We find no strong evidence to support an interaction between NGC 4649 and the neigh¬ 
bouring spiral galaxy NGC 4647. We identify interesting GC kinematic features in our data, 
such as counter-rotating subgroups and bumpy kinematic profiles, which encode more clues 
about the formation history of NGC 4649. 


1 INTRODUCTION 

Galaxies are the building blocks of the visible Universe and the 
best tools to study its structure. Galaxies are alike in many ways, 
suggesting that similar underlying mechanisms regulate their for¬ 


mation, but they are not identical to each other, impl ying that for¬ 
matio n processes vary slightly from galaxy to galaxy dUintott et al.l 
l2008h . While high redshift galaxy surveys can be used to infer “av¬ 
erage” formation mechanisms for a sample of galaxies, studying 
nearby galaxies in detail can allow us to reconstruct their particular 
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and unique formation histories. Moreover, observations of galaxies 
at 2 = 0 provide t he end product that computer simulations attempt 
to rep roduce te.g.. lHoffnian et al.l2010l : IWu et alJ2014l : lNaab et al.| 


to rep rc 
l2014h . 


In th e standard galaxy formation theory, small structures col¬ 
lapse first dWhite & Reeslll97^ : ISearle & Zinnll97^ : IZolotov et al.l 
l2009l) . and then grow hierarchically into larger structures via 
galaxy mergers and / or via accretion of sate l lite galaxies (e.g.. 


Logez^Sanjuanrt^ 20 id : iFont et aP 1201 ll : Ivan Dokkum et all 
2014 Tasca et al.ll20i% The latter is responsible for the building- 


up of galaxy outskirts (which we will refer to as stellar haloes) from 
high (2 « 2) to low redshift l lTal & van Dokku^l201 lUOser et al] 
l20ld : iKhochfar et al.ll201 ill . The various processes which shaped 
galaxies with time can be simulated and compared with obser¬ 
vations of galaxies at a particular redshift. For example, minor 
mergers can c reate stellar shells and tidal s treams observable in 
deep imaging jTal et aD l2009l: lEb rov 3 I 2 OI 3 I : lAtkinson et ahll^oT^ : 
IDuc et ^ 2oiar whereas the secular accretion of satellite galax- 
ies can be studied with photometry of stacked galaxies (e .g., 

lOser et aI]|20ld : lTal & van DokkunJl201 ihlP’Souza et alj20l4 or 

by studying the chemistry and kinematics of stellar haloes (e.g .. 


iForbes et'^l201 ll : iRomanowskv et alj2012l : ICoccato et alj|2013 

Galaxy haloes are difficult to study because they are opti¬ 
cally faint. Globular clusters (GCs), on the other hand, are much 
more observationally convenient for studying galaxy haloes. They 
are relatively easy to detect because of their high surface bright- 
ness and they generall y populate haloes in large numbers (e.g., 
ISrodie & Stradeill2006ll . The outermost GCs can be used to map 
out halo properties (e. g., kinematics and metallicity) out to tens 
of effective radii ( e.g..]ostroy et 1993: Schuberth et alj|20l^ : 
lusher et al.|[20ll: IPota et alJl2013l : yForbes et al.ll201 ih . Their old 
ages and their direct connection with the star forming ep isodes 
in a galaxy’s history jStrader et alj|2003 : IPuzia et alj|2005t) mean 
that they are ideal t racers of the assembly processes that shaped 
the host galaxy (e .g., Romanowsk^et_^ 2012l : lLeaman et al.l2013l : 


iFoster et alj2014l:lvelianoski et al.ll2oT4). 

In this paper we study the GC system of the E2 early-type 
galaxy NGC 4649 (M60), the third brightest galaxy in the Virgo 
cluster. NGC 4649 is itself at the centre of a small group of g alax¬ 
ies an d is the dominant member of the galaxy-pair Arp 116 
1 19661) . The proximity to the disturbed spiral galaxy NGC 4647 (13 
kpc in projection) potentially makes NGC 4649 an example of a 
pre major-merger between two massive galaxi es in the local Uni¬ 
verse, although this connec tion is still debated l lYoung et al.ll2006l : 
Ide Griis & Robertsonll2003) . 

Galaxies like NGC 4649 are ideal laboratories to test galaxy 
formation models. This galaxy has been scrutinized from dif¬ 
ferent angles, revealing the portrait of a prototypical red and 
dead elliptical galaxy w ith a central super massive black hole 
jShen^_Gel^rdj l2010t) . dark matter halo jBridges et alj l20Q6l: 
DasetalJ|201T ) and X-ray halo 1 O’Sullivan et alj 2003 ; Das et al.l 


201(]l:lHumphrev et alj2008l:IPaggi et alj2014l). 

Recent space-based observations of NGC 4649 JStrader et al.l 
I 2 OI 2 I: ILuo et al.ll20l3l: iNorris et al.ll2014l) . suggest that this galaxy 
interacted with a smaller galaxy with mass M ~ 10^°Mq 
JS eth et alj|2014l) . This event should have deposited 100 GCs and 
many m ore stars into the h alo of NGC 4649 during the course of 
its orbit jHarris et al.ll2013l) . Indeed, asymmetries have been found 
in the two-dimensional distri butions of GCs and low-mass X-ray 
binar ies around this galaxy jMineo et al. I I 2 OI 4 I : iD’Abrusco et al.l 
l2014h . but these have not yet been linked to any particular inter¬ 
action that occurred in NGC 4649. 


In addition to recent space-based observations, the GC sys- 
tem of NGC 4649 has been studied with ground-based telescopes 

(|Bridges_et al.l2006l:IPierce et al.l2006l : lLee et al.l2008l : lFaifer et al.l 

I 2 OIII : Chies-Santos et al. 201^. A spectroscopic follow up of 


121 GCs I Hwang et al. 20081) revealed that the GC system of 


NGC 4649 has a large rotation amplitude. This is a rare feature 


sure supported with negligible or weak rotation (e.g.. Cote et al. 

2003; 

Bergond et al.l l2006l: ISchuberth et al.l 

l2010bl; Strader et al. 

2011; 

Norris et alj 20121; IPota et alj 2013; 

Richtler et al. 12014). 

although exceptions exist (e.g.. Puzia et al 

2004 Arnold et al. 

201ll;ISchuberth et al.ll2012l:lBlom et alJl20U 

). On the other hand. 

Bridges et alJ (2006h obtained radial velocities for 38 GCs in 


NGC 4649 and found no rotation, probably because of their rel¬ 
atively small sample size. 

In this paper we present the largest spectro-photometric cata¬ 
logue of GCs around NGC 4649. We construct a wide-held pho¬ 
tometric GC catalogue by combining literature (HST), archival 
(CFHT) and new (Subaru) observations. We follow-up with spec¬ 
troscopy of hundreds of GC candidates using joint observations 
from three multi-object spectrographs mounted on the Keck, MMT 
and Gemini telescopes. Our new spectroscopic ca talogue is a fac- 
tor of 3.5 larger than the current literature dataset dUee et al.ll2008t) 
and has a factor of 3 greater velocity accuracy. This study exploits 
datasets from two complementary sur veys of extragalactic GCs: the 
SLUGGS survey dBrodie et al.ll2014l) . and an ongoing survey car¬ 
ried out with Gemini/GMOS (Bridges et ah, in preparation). These 
two surveys combined, contribute 90 per cent of the hnal sample 
of conhrmed GCs. We use the photometric and kinematic proper¬ 
ties of the NGC 4649 GC system (including the blue and red GC 
subpopulations) to study the formation history of NGC 4649. A 
follow-up paper (Gebhardt et al., in preparation) will make use of 
the dataset presented in this work to model the mass content of 


NGC 4649. __ 

We adopt a distance of 16.5 Mpc dMei et al.ll2007l) . an effec¬ 
tive radius Re — 66 arcsec = 5.3 k pc, an axis ratio q = 0.84 
and a position angle PA — 93 deg dBrodie et alj|2014l) . Galac- 
tocentric distances R are expressed though the circularized radius 
R = ^ X'^q + jq, where X and Y are the Cartesian coordi¬ 
nates of an object with NGC 4649 at the origin. The absolute mag¬ 
nitude of NGC 4649 is Mb = -21.47 mag or Mv = -22.38 
mag. 


This paper is structured as follows. Imaging observations and 
analysis are presented in Sections |^ and respectively. The pho¬ 
tometric results, including colour gradients and GC surface den¬ 
sity are presented in Section|4] The spectroscopic observations and 
their outcomes are discussed in Section and Section (6] respec¬ 
tively. We next explain the steps needed to build up the spectro¬ 
scopic GC master catalogue: repeated GC measurements (Section 
|7]l and colour and velocity uncertainty calibration (Sectionj^. We 
give an overview of the spectroscopic sample in Section and ex¬ 
plain the tools to quantify GC kinematics in Section[T0] The kine¬ 
matic modelling results are given in Section EU We discuss our 
results in Sectionfl^and summarize the paper in SectionfOl 


2 IMAGING OBSERVATIONS 

Imaging is the hrst step in identifying extragalactic GCs. It allows 
us to identify and characterize the photometric properties of the 
NGC 4649 GC system, and also to prune out the bulk of the con- 
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The globular cluster system ofNGC 4649 3 



Figure 1. Ovei'view of the observations. Shown is a Digitized Sky Survey image ofNGC 4649 and suiToundings. The field of view is 1 deg^, which corresponds 
to 290 kpc X 290 kpc at this distance. NGC 4649 is the elliptical galaxy at the centre of the image, with the spiral NGC 4647 visible on the immediate right. 
Other galaxies in the field are labelled. Black dashed polygons represent the fields-of-view of the imaging cameras (labelled accordingly). Coloured solid 
outlines enclose the regions in the sky mapped by our three multi-object spectrographs (labelled). A 10 arcmin~ 48 kpc scale-bar is shown at the bottom of 
the image. North is up. East is left. 


taminants (Galactic stars and background galaxies). Imaging is also 
necessary for selecting objects to be followed up spectroscopically. 

Our imaging comes from three sources: HST/ACS, 
CFHT/MegaCam and Subaru/Suprime-Cam. The fields of view of 
the three instruments, along with the footprints of the multi-object 
spectrographs, are shown in Figure[T] CFHT/MegaCam covers the 
field of interest, and it was supplemented with Subaru/Suprime- 
Cam data to gain extra coverage on the East of NGC 4649 and to 
test for systematics in our ground-based photometry. 

2.1 HST/ACS 

The HST imaging is discussed in istrader et alj ( 1201 2h and we refer 
the reader to this paper for a detailed description of the data. The 


dataset consists of six ACS pointings, which extend up to 6 arcmin 
from the galaxy center. The layout of the pointings is such that it 
avoids the foreground spiral NGC 4647. GC selection is based on 
{g — z) colours, z magnitudes and half-light radii r^. The latte r 
can be measured down to ~ 0.1 x FWHM dSnitler et al.ll2006h . 
which corresponds to ~ 0.7 pc at the distance of NGC 4649. Unre¬ 
solved sources were not selected as GCs. The latter is the physical 
radius which contains half of the total light of the object. The final 
HST/AC5 catalogue consists of 1603 GC candidates. 

2.2 CFHT/MegaCam 

The area surrounding NGC 4649 has been imaged with 
CFHT/MegaCam as part of the Next Generation Virgo Survey 
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dFerrarese et alj 1^0121) . MegaCam has a field-of-view of 0.96 x 
0.94 deg^ and a pixel scale of 0.187 arcsec/pixel. Reduced 
CFHT/MegaCam images and the respective weight-images in 
gri filte rs were down loaded from the CFHT/MegaPipe legacy 
archive dOwynlllOOSh . Images are provided already calibrated to 
standard SDSS filters. The logs of the CFHT observations are given 
in Table[T] The seeing for all images is < 0.8 arcsec. 

We use a total of five MegaCam pointings. Four of these are in 
Mgiz-filters, but they are not centered on NGC 4649. A fifth point¬ 
ing in the r-filter is centered on NGC 4649. The final CFHT field 
of interest spans roughly 1 degree in declination and half a degree 
in right ascension, as shown in Figure [T] We discard 30 arcmin in 
right ascension as we are not interested in the sky to the West of 
NGC 4621. 

2.3 Subaru/Suprime-Cam 

NGC 4649 was imaged with Subaru/Suprime-Cam in gri filters. 
The field-of-view of Supri me-Cam is 34 x 27 ar cmin^, with a pixel 
scale of 0.2 arcsec/pixel dMivazaki et alj|^002h . The camera was 
pointed at roughly the centre of NGC 4649. A summary of the Sub¬ 
aru observations is given in Table [T] The seeing was sub-arcsec for 
the g and i images, and 1.1 arcsec for the r image. 

We use a custom made pipeline for reducing the Sub- 
aru Suprime-Cam imaging, based on the SFRED2 pipeline of 
lOuchi et aP d2004h . The reduction consists of the following steps: 

• We subtract the bias measured from the overscan regions for 
each chip. We create the flat field frames by me dian combining the 
dome-flats for each chip. We use L.A. Cosmic d vanDokkumll200 11) 
to remove the cosmic rays from the object frames. 

• We correct for chip-to-chip quantum efficiency variations and 
scale individual chips for any time-variable gain. We correct for 
the atmospheric distortion, mask out the auto-guider probe in each 
frame and sky subtract using a simple median of medians. 

• We use a Monte Carlo sampling method to efficiently find the as¬ 
trometric solution which minimizes the positional offsets between 
SDSS and point sources in our data. 

• We multiplicatively scale each exposure to the median through¬ 
put level of all the exposures before co-adding them into a final 
image. The final mosaic image for each filter was performed with 
the software Montage Q. 


3 IMAGING ANALYSIS 

3.1 Catalogue extraction and calibration 

We use SExtractor dSertin & Amoutsll996h to extract sources from 
the Suprime-Cam and MegaCam images. We extract sources from 
the two cameras separately, and we combine them later on. We run 
SExtractor on all the fields listed in Table Q] Only sources above 
2.5 (t are extracted. The zero point magni tude is set to 30 mag for 
the CFHT images dFerrarese et alj|2012ll . and to 25 mag for the 
Subaru images. The weight-maps downloaded from the MegaPipe 
archive were used to optimise the extraction of MegaCam objects. 
The weight-map for the Suprime-Cam dataset was set to “back¬ 
ground”, meaning that we let SExtractor compute the variance map 
from the science images themselves. 

Eor each extracted source, we measure: a) coordinates; 

^ http://montage.ipac.caltech.edu/index.html 


Obs ID-Filter 

Date 

Exposure time 

Seeing 


[HST] 

[seconds] 

[arcsec] 

Subaru-^ 

2010-04-11 

1160 

0.8 

Subaru-r 

2010-04-11 

1260 

1.1 

Subaru-i 

2010-04-11 

490 

0.7 

G008.190.711+11.596-r 

2008-03-01 

2160 

0.8 

NGVS-t3-t0-g 

2009-02-27 

3170 

0.6 

NGVS-t3-t0-i 

2009-02-27 

2055 

0.6 

NGVS-t3-l-g 

2009-02-27 

3170 

0.6 

NGVS-t3-l-j 

2009-02-27 

2055 

0.6 

NGVS-t2-tO-g 

2009-02-25 

3170 

0.6 

NGVS-t2-tO-/ 

2009-02-15 

2055 

0.6 

NGVS-t2-l-g 

2009-05-25 

3170 

0.6 

NGVS-t2-l-i 

2009-06-20 

2055 

0.6 


Table 1. Ground-based photometric observations. Listed are the telescope 
and filter, the observation date in Hawaiian Standard Time (HST), the ex¬ 
posure time and the seeing. The overall seeing conditions are excellent. 


b) magbest, which we use as our primary magnitude measure¬ 
ment; c) a set of magnitudes computed within different apertures 
(from 4 to 11 pixels in diameter); d) structural parameters, such as 
the ellipticity, position angle and EWHM. 

The photometric properties of the objects in common between 
the different g and i MegaCam pointings are in good agreement 
with each other. Selecting objects brighter than i = 23 mag, we 
find the median magnitude difference to be Ag = 0.006 mag and 
Ai = 0.008 mag, respectively. The extracted properties for these 
objects are averaged together so that every object is unique, with 
magnitude measurements in all three filters. 

Point sources are selected based on the difference between the 
magnitude measured within 4 pix and within 8 pix. This difference 
is small for point sources (typically < 0.6 mag) because their light 
profiles fall off more rapidly than those of extended sources. The 
latter become important source of contamination for {i > 22 mag). 

With the point-source catalogues in hand, we calibrate the 
magnitudes to standard SDSS gri filters. We select a bright sub¬ 
sample (18 < i < 21 mag) and match these sources to the SDSS 
DR7 point-source catalogue, which covers the entire field-of-view 
in Pigure[T] 

The MegaCam dataset was already calibrated to SDSS mag¬ 
nitudes and, in fact, we find the zero-point corrections to be al¬ 
ways < 0.05 mag, within the photometric uncertainties. Eor the 
Suprime-Cam dataset we find: grsoss — 5 sub = 3.46 ± 0.05 mag, 
fSDSS — fsuB = 3.76 ± 0.05 mag, /sdss — *sub = 3.66 ± 0.04 
mag, respectively. 

3.2 Comparison between CFHT and Subaru 

Before merging the MegaCam and Suprime-Cam datasets, we 
study how the magnitudes extracted from both cameras compare 
with each other. We match the two catalogues selecting only bright 
(19 < i < 23) objects outside R = 300 arcsec, because ground- 
based imaging is notoriously incomplete near the centres of the 
galaxies. 

Eigurel^shows the magnitude difference in three bands of the 
~ 1800 point sources in common between MegaCam and Suprime- 
Cam. We are mostly interested in objects with 0.65 {g — i) ^ 1.4 
because this is the colour range populated by extragalactic GCs, as 
we will discuss below. In this specific colour-magnitude range, the 
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Figure 2. Comparison between Subaru and CFHT datasets. Top, central and 
bottom panel shows the magnitude difference in g, r, i filters, respectively, 
of the point-sources in common between Subaru and CFHT. The green line 
marks zero magnitude difference. Only objects with 19 < i < 23 and 
0 < (g — i) < 2.5 are shown. Globular clusters fall in the colour range 
0.07 ^ {g — i) ^ 1.4, which is highlighted with black points. In this colour 
range, the CFHT and Subaru datasets agree with each other. The deviation 
of ~ 0.07 mag for g and i filters, and ~ 0.09 mag for the r band. 


root-mean-square (rms) of the magnitude difference is rms(Ag) = 
0.070 mag, rms(Ar) = 0.095 mag, rms(Ai) = 0.064 mag for 
the g, r, i filters respectively. The larger value of rms(Ar) is due 
to the relatively poor quality of the Suprime-Cam r-band image. 
No significant trend with {g — i) colours is observed in the colour 
range under investigation. We conclude that the magnitudes mea¬ 
sured from MegaCam and from Suprime-Cam are in good agree¬ 
ment. 

We then combine the point-source catalogues from Mega- 
Cam and Suprime-Cam. If an object is present in both datasets, 
we weight-average its magnitudes in the three filters. If an ob¬ 
ject is present in only one catalogue, it is added to the final mas¬ 
ter catalogue without applying any magnitude correction. Lastly, 
we correct the final gri magnitudes for Galactic extinction using 
Ag = 0.087 mag, Ar — 0.0 60 mag, Ai = 0.045 mag, respec¬ 
tively JSchlaflv & Finkbeineill201 Ih .These values vary by 10 per 
cent across our field-of-view of interest. The combined MegaCam- 
Suprime-Cam point-source catalogue consists of ~ 20500 unique 
objects, including interlopers. 


3.3 Globular cluster selection 


We use two steps to extract bona-fide GCs from the point-source 
catalogue. We first exploit the HSTlACS catalogue to check where 
GCs lie in colour-colour and colour-magnitude space, and then we 
select GCs based on this comparison. 

GCs are partially resolved in HST i mages. In NGC 46 49 the 
average GC half-light radii are rj, « 2 pc jStrader et al.l2Q12h. sim - 
ilar to what is found in other galaxies (e.g.. Masters et al.ll201oh . 
This means that even without spectroscopic confirmation, we can 



0.4 0.6 0.8 1.0 1.2 1.4 1.6 


(g - i) [mag] 


Figure 3. Colour-magnitude and colour-colour magnitude diagrams. For 
clarity, only ground-based sources within 600 arcsec from NGC 4649 are 
shown (small grey points). Large black points are partially resolved GCs 
from H ST. Orange open squares are the confirmed GCs from iHwang et alj 
I 2 OO 8 I) . The bottom panel shows the colour-colour magnitude plot. The cen¬ 
tral panel is the colour-magnitude plot. Horizontal lines indicate the TOM 
of the GC luminosity function i (TOM) = 23.1 mag, and the magnitude 
of oj Cen, i = 20 mag. Objects inside the black lines are flagged as GC 
candidates. The top panel is the numerical histogram of the selected GCs, 
with the best double-Gaussian fit to the data. 


flag objects with apparent rn < 10 pc as GCs likely associated with 
NGC 4649, although some contamination is expected at faint mag¬ 
nitudes. From the ACS catalogue, we select GC candidates brighter 
than 2 = 23 mag and with R > 100 arcsec and we match them 
with our combined ground-based catalogue. 

We find roughly 500 objects in common. Then we look at the 
colour-magnitude and colour-colour diagrams of the point sources 
with ground-based photometry, and we highlight these 500 GCs 
with size confirmation. The result is shown in Figure We also 
matched the AC S catalogue with the 121 spectroscopically con¬ 
firmed GCs from iHwang et al.l ( 1200 Sl) . which will be discussed in 
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Section|5] We found 84 objects in common, which are also shown 
in Figure [3 

NGC 4649 GCs populate a well defined area of both the 
colour-magnitude and colour-colour diagrams. In the central panel 
of Figure we see the well known GC colour bimodality and the 
blue tilt, which causes t he most luminous blue GCs to bend to¬ 
wards redder colors te.g.- ISnitler et al.l2006l : lHmis et al .l2006h . At 
the faint and bright GC selection limits we show, respectively, the 
expected turn-over magnitude ( TOM) of the GC lu minosity func¬ 
tion (Mtom,* = —7.97 mag: Llordan et al.ll2007h and the mag¬ 
nitude of uj Cen, the brightest star cluster in the Milky Way, for 
which we assu me a magnitude Mj = — 11.0 mag (i = 20.1 mag) 
obtained from Vanderbeke et alj ( l2014h adopting the distance of 
Ivan de Ven et alJ(l2006h . A handful objects are indeed brighter than 
this magnitude. We will show that some of these are ultra compact 
dwarfs (UCDs) with spectroscopic confirmation (see 

On the bottom panel of Figurel^we notice the same grouping 
of GCs in a narrow colour-colour range. The tail of objects at very 
red colours are Galactic red giant stars which enter the GC main 
sequence at {g — i) ~ 1.4 mag and contaminate the sample. The 
objects at {g — i) ^ 0.6 are Galactic blue young stars. Roughly 20 
per cent of the GC candidates lie outside the selection box. This is 
due to photometric uncertainties and to the poor quality of ground- 
based photometry near the centre. Selecting only objects with z < 
22 mag and with R > 200 arcsec, the fraction of outliers reduces 
to 3% . 


Objects falling in the colour range highlighted in Figurel^are 
flagged as GC candidates. Objects falling outside the selection box 
by less than la are also flagged as GCs. We require a GC to be 
brighter than the turn-over-magnitude (i < 23.1 mag) and fainter 
than (jj Cen. These selection criteria returned ~ 4000 GC candi¬ 
dates, with a 14 per cent contamination rate based on the back¬ 
ground value which will be computed in Section|4] 

The colour distribution of the selected GCs is bimodal (top 
panel of Figure^- To quantity this bimodality, we use th e Kaye’s 
Mixture Model (KMM) algorithm l lAshman et al.l 1 199^ to fit a 
double-Gaussian to the GC colour distribution. This calculation 
was performed only on objects within 600 arcsec from the galaxy 
centre in order to prune the bulk of contaminants which populate 
the outer regions of the galaxy (see Section|4ll. 

We find bimodality to be statistically significant with a p-value 
< 10~^. The blue GC subpopulation is fitted by a Gaussian which 
peaks at p(p — i) = 0.83 mag, with a dispersion of (T(p — i) = 0.07 
mag. For the red GC subpopulation, we find fj,(g — i) = 1.08 
mag , with a dispersion of a{g — i) = 0.12 mag. For compari¬ 
son, |Faiferer^(l2oTTl) found the following values based on Gem- 
ini/GMOS photometry: fi{g — i) = 0.78, a{g — i) — 0.09, and 
P(5 ~ = 1.08, a{g — i) = 0.11, for the blue and red subpopu¬ 

lations, respectively. The local minimum of the combined double- 
Gaussian is at {g — i) = 1.0 mag. This is the value adopted to 
separate blue and red GCs in our ground-based imaging. 


4 PHOTOMETRIC RESULTS 

We now use our photometrically-selected GCs to study two impor¬ 
tant relations: how the GC number density and the GC colours vary 
with radius. The first is an essential observable for dynamical mod¬ 
elling of galaxies, but it is also needed to quantify contamination 
from interlopers. The shape of the second relation is an important 
prediction of theories of the hierarchical growth of galaxies. 


R[kpc] 



Figure 4. Globular cluster surface density profiles. Blue and red GCs are 
colour-coded accordingly. Open and filled circles are from HST/ACS and 
from ground-based photometry, respectively. The two solid lines are the 
best fits to equation[T] Green squares are PNe data. The black line is the ar¬ 
bitrarily re-n ormalized stellar surface brightness profile of NGC 4649 in the 
U-band from iKormendv et al.l i2009h . The magenta thick and dashed lines 
are the outer slope of the X-ray surface brightness profile from ROSAT and 
from Chandra, respectively. The red GCs are more centrally concentrated 
than the blue GCs. 


4.1 Globular cluster surface density 

We divide our sample into blue and red GCs at {g — i) = 1.0 
mag. We count the number of GCs in elliptical annuli centered on 
NGC 4649 and we divide this number by the area of the annulus. 
Saturated stars, surrounding galaxies and chip gaps in the CHFT 
images were masked out. For each bin, we compute the Poissonian 
error using y/N / Area, where N is the number of objects per bin. 

As our ground-based imaging is incomplete near the galaxy 
centre, we s upplement our data with the ACS GC surface density 
profile from iMineo et alj ( l2014h (their Figure 2), which covers the 
innermost 5 arcmin of the galaxy. This was obtained selecting GCs 
brighter than 2 < 22.2 mag and it was renormalized to match the 
galaxy surface brightness. Therefore, the ACS surface density pro¬ 
file does not represent absolute counts of objects with 2 < 22.2. We 
construct the ground based surface density profile selecting GCs 
about 1 magnitude fainter with respect to ACS. To account for this 
normalization difference, we count the GCs selected with our cri¬ 
teria in the radial range 150 < R < 270 arcsec , and we compare it 
to the number of GCs from iMineo et al.l ( l2014h (their Figure 2) in 
the same radial range. We find that a factor of 4.3 and 5.7 is needed 
to match the ACS-based surface density profile for the blue and red 
GCs respectively, with our ground-based surface density profiles. 

The combined GC surface density profiles for the blue and red 
GCs are shown in Figure|4] 

At first glance, our results show the typical spatial features 
seen in other GC systems: the red GCs are more cen trally concen¬ 
trated , while the blue GCs are more extended (e.g.. lBassino et all 
1 20061) . The flattening of the number density at large radii is due 
to contaminants (faint Galactic stars and unresolved galaxies), and 
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marks the detectable boundary of the GC system. To quantify the 
background level hg, assumed to be constant across the entire im - 
age, we fit the blue and red GC data with a modified lSersi 3C96i: 


Ngc{R) 


X exp 




+ bg, ( 1 ) 


where n is the Sersic index and Ne. is the surface density at the 
effective radius R^. The parameter is linked to n via = 
1.9992n - 0.3271. 

The best fits to equation [T] are shown in Figure |4] For the red 
GCs, we find Rf, = 4.1 ± 1.1 arcmin = 20 ± 7 kpc, n = 2.5 ± 0.5 
and bg — 0.2 ± 0.1 objects/arcmin^. For the blue GCs, we find 
Re = 12.6 ±7.1 arcmin = 60 ± 34 kpc, n = 2.4 ± 0.6 and 
bg = 0.23 ± 0.19 objects/arcmin^. Therefore the blue GC system 
is a factor of three more extended than the red GCs. 

In Figure [4] we also compare the GC spatial distribution 
with the surface density of the following tracers; (i) the stel- 
lar surface brigh t ness p rofile of NGC 4649 in the V-band from 
iKormendv et al.l ( 120091): (ii) the X-ray su rface brightness pro¬ 
file from ROSAT lo Sullivan et alj 1200 3t) and from Chandra 
dPaggi et alj|2014h . In this case, we show a /3-model best fit to 
X-ray data. The best fit is shown up to the radius where the X- 
ray emission drops below the background level (~ 4.8 arcmin) for 
the ROSAT data, (iii) Planetary nebulae (PNe) surface density data 
from lPas et alj ) |201 ih . with the caveat that these were obtained se¬ 
lecting objects within a cone aligned with the galaxy major axis. 
The three tracers were arbitrarily re-normalized in order to be com¬ 
pared with the GC surface density profiles. 

We can see that the slope of the stellar and X-ray surface 
brightness is qualitatively similar to the slope of the red GCs, with 
the caveat that this similarity is dependent on the arbitrary nor¬ 
malization. X-r ay haloes are usually more extended than the stel- 
lar component dSarazin et al.ll200ll ; lForte et al.ll2003 : lForbes et al.l 
l2012h . Results from Chandra and from ROSAT are consistent with 
each other. The similarity between stars-PNe and red GCs is ex¬ 
pected if the galaxy bulge an d the red GCs form ed in a similar 
fashion at similar epochs (e.g.,[Sha 2 iro_etalJ 201 ^. 

We follow the approach of Forbes et alj d2012h to quantify the 
similarities between the distribution of GCs and X-ray gas. We fit 
the GC density profiles in Figure |4] with a /3-model, which is a 
cored-power law described by 7(r) oc where /3 is the 

power-law slope we are interested in. The background values bg 
derived above were taken into account d uring the fit. A fit to th e 
ROSAT data returns Px = 0.56 ± 0 .08 dO’Sullivan et al.ll200^ . 
Using new deep Chandra data from IPaggi et alj d2014l) . one ob¬ 
tains consistent results with Px = 0.52±0.01. For the red and blue 
GCs, we find pn = 0.48 ± 0.03 and Pb = 0.44 ± 0.05, respec¬ 
tively. In this case the larger uncertainty reflects the uncertainty on 
the GC background value. This exercise shows the outer slope of 
blue and red GCs are consistent within the errors, but the red GC 
profile is slightly steeper and more similar to the X-ray slope. 



0.7 0.9 1.1 1.3 
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Figure 5. Photometric calibration between ground-based and space-based 
observations. Data points represent partially resolved GCs in HST images, 
with constraints on galactocentric distance and magnitude (see text). The 
red line is a weighted least-squares fit to the data. 


R/Re 



R [kpc] 

Figure 6. Radial variation of GC colours. Grey and black small points are 
GC candidates from HST observations and from our ground based imaging, 
respectively. Black lines are best fits to the data up to 60 kpc (see text). 
Blue and red filled large circles show the colour peaks from KMM for a 
given radial range. Blue and red open large circles represent the peaks of 
the Gaussian-smoothed colour histogram in each bin (see text). The green 
line is the stellar metallicity profile converted to (g — z) (see text). The GC 
subpopulation colour gradients ai‘e steeper towards the galaxy centre. 


4.2 Globular cluster colour gradients 

We study the radial variations of GC colours by combining HST 
and ground-based imaging. It is more efficient to calibrate both 
datasets to a common photometric colour and then merge them. 
We choose to transform the ground-based [g — i) to ACS {g — z) 
because this can be then be co nverted to metallicity through the 
conversion o f lPeng et Tl i2006l) . 

As in the previous section, we match the objects in common 


between HST and the combined Subaru±CFHT catalogue. We re¬ 
quire ACS sources to be partially resolved with sizes < 10 pc 
and magnitudes 19 < i < 22 mag. We only consider sources with 
galactocentric radius 100 < R < 400 arcsec because the ground- 
based imaging is unreliable within 100 arcsec. 

We find 241 objects obeying the above criteria. A weighted 
last-squares fit to the data returns the colour transformation we are 
interested in: 

(g-z) = (-0.18 ± 0.03) ± (1.36 ± 0.03)(5 - i) (2) 
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The best fit, along with the data, is shown in Figure [S] The rms of 
the data with respect to the best fit line is 0.08 mag. The quality 
of the fit is satisfactory. We also convert the ACS g magnitudes to 
ground based g magnitudes. We find (?acs = (0.62 ± 0.20) + 
(0.97 ± 0.01)(; with an rms scatter of 0.09 mag. Similarly, we find 
^Acs = (0.07 ± 0.25) + (0.99 ±0.01)zACS, with a scatter of 0.13 
mag. We apply eq. to our ground-based GC catalogue and we 
merge this into the ACS catalogue. Repeated objects are counted 
once and priority is given to objects with ACS photometry. 

We construct the colour gradient s by selecting (21 < 2 < 23) 
(consistently with lStrader et alj2013h to avoid contamination from 
interlopers at faint magnitudes and from the blue-tilt at bright mag¬ 
nitudes. We also clip objects within 5 arcmin of large galaxies in 
the field to avoid mixing their GC systems with that of NGC 4649. 
The results are shown in Figure |6] 

We can see that the ACS and ground-based datasets agree and 
that they are complementary to each other. We study the colour pro¬ 
files out to ~ 16 arcmin ~ 80 kpc, as dictated by the background 
estimate from equation[T] This radius corresponds to 15 Re, where 
Re = 66 arcsec is the effective radius of the diffuse starlight. 

We compute GC colour radial profiles using two different 
methods. First, we bin the data radially with a 5 kpc bin size for 
the innermost 20 kpc, and with a 10 kpc bin size for R > 20 
kpc. For each bin, we fit a double Gaussian with KMM, and de¬ 
rive the blue and red colour peaks of the fitted Gaussians. Second, 
we smooth the colour distribution with a Gaussian kernel of 0.05 
mag and we find the blue and red peaks of the resulting smooth 
distributions. This second method is model independent, whereas 
the first method requires assumptions that may not provide a re¬ 
alistic representation of the data. The results are shown in Figure 
as filled points for the first method (KMM) and open points for 
the second method (smooth kernel), respectively. We also compare 
our re sults with the stellar metallicity profile from IPastorello et aP 
l 2014h. converted from [Z/H] to {g — z) using equation B1 from 
lUsher et^ l l2012l) . 

We can see that the results in Figure|6]are not very sensitive to 
the method adopted to compute the colour peaks. Both the blue and 
the red GCs show negative g radients within « 20 kpc, as already 
found in IStrader et al.l j2012h . Outside this radius, the colour gra¬ 
dients are consistent with being constant. The overall slope of the 
stellar colour profile is in good agreement with the results for the 
red GCs. This provides further evidence that stars from the galaxy 
spheroid and red GCs were assembled in a similar fashion, at least 
in the very central regions. 

We quantify the gradients by fitting the function {g — 
z) = alog{R/Re) + b, where Re = 5.3 kpc is the stellar ef¬ 
fective radius of NGC 4649. We fit all individual GC datapoints in 
Figurel^up to 60 kpc to minimize contamination from background 
sources. We divide the two GC subpopulations by adopting a cut 
at (g — z) > 1.2 for the red GCs and {g — z) < 1.1 for the blue 
GCs, respectively. This is done to avoid mixing between the two 
GC subpopulations at intermediate colours. For the red GCs, we 
find a = -0.07 ± 0.01 and b = 1.44 ± 0.01. For the blue GCs, 
we find a = —0.005 ± 0.010 and b = 0.94 ± 0.01, indicating an 
overall shallower colour gradient with respect to the red GCs. 

The red GC colour gradient is wiggly in the range between 
10 and 20 kpc, suggesting the presence of s ubstructures in the data 
jStrader et al.l2012 : lD’Abrusco et al.l2014l) . and it is very steep in¬ 
side 10 kpc. This is due to a group of GCs with very red colours 
{g — z) > 1.5, which populate only the innermost 10 kpc of the 
galaxy. 


5 SPECTROSCOPY 

Our spectroscopic sample comes from three multi-object spectro¬ 
graphs: Keck/DEIMOS, MMT/Hectospec and Gemini/GMOS. The 
footprint of the observations of each instrument is shown in Fig¬ 
ure [T] The spectral coverage of the three spectrographs is diverse, 
but some spectral regions overlap. Figure [7] shows one of the nine 
GC spectra observed by all three spectrographs. In this section we 
describe the acquisition and reduction of the spectroscopic data. 


5.1 Keck/DEIMOS 

DEIMOS is a mul ti-object spectrogr aph mounted on the 10 meter 
Keck II telescope l lFaber et al.ll2003h . One DEIMOS mask covers 
15x6 arcmin^ on the sky and it can contain up to 150 slits. 

A total of four DEIMOS masks were observed around 
NGC 4649. Three masks were observed on 2012 Jan 17, and the 
fourth mask was observed on 2013 Jan II. Slitmasks were designed 
giving priority to ACS GC candidates with size measurements. If 
no space-based photometry was available, slits were assigned to 
bright sample of GCs selected from the CFHT imaging discussed 
in Section Suprime-Cam images were not yet available at the 
time of the mask design. 

All observations were performed in sub-arcsec seeing condi¬ 
tions. The exposure time was 2 hours per mask. The instrument 
setup consisted of a 1200 1/mm grating, a dispersion of 0.33 Aj 
pixel, a slit width of 1 arcsec and the OG550 filter. With this setup 
the GC spectra cover the region from the Ha line (6563 A) to the 
Calcium Triplet (CaT) feature (8500-8700 A), with a median ve¬ 
locity accuracy of ~ 15 km s~^. The spectra were reduced with 
a dedicated pipeline, which returns calibrated, sky-subtracted one¬ 
dimensional spectra dCooner et al.ll20l3 : lNewman et al.ll2013h . 

GC radial velocities were measured with iraf/fxcor using thir¬ 
teen stars of different spectral type as templates observed using the 
same set-up. The average velocity scatter between different stellar 
templates is ~ 3 km s“^ for all four DEIMOS masks. The Doppler 
shift was computed solely from of the CaT lines, because the Ha 
line was not always covered with our instrument setup. The final 
GC radial velocity was taken as the average of the values from the 
thirteen stellar templates. The redshift-corrected spectra were in¬ 
spected by eye to check that the CaT lines laid at the rest frame 
wavelength. The error budget on the velocity takes into account 
both the internal error given by fxcor and the scatter due to differ¬ 
ent stellar templates. 


5.2 Gemini/GMOS 

GMOS is a multi-object spec trograph and imag er mounted on the 
8 -meter GEMINI telescopes faook et alj|2004h . The GMOS cam¬ 
era has three 2048 x 2048 pixel CCDs with a scale of 0.07 arc- 
sec/pixel and a field of view 5.5 x 5.5 arcmin^. GMOS was used to 
obtain spectra for globular clusters in five fields in NG C 4649 (see 
Figurj T). GC selection was performed as described in lEaifer et al.l 
( l201lh . GC spectroscopy for a field positioned near the galaxy cen¬ 
ter (Field 1) was obtained in May 2003 on Gemini North a nd was 
presented in [Bridges et all j200^ and IPierce et aP (l2006h . while 
spectroscopy for four additional fields (Fields 2-5) was obtained 
in the spring of 2010 on Gemini South and presented in this paper. 

GC candidates for follow-up spectroscopy were selected using 
Sloan g' and i' images taken with GMOS-N in April 2002. Images 
were combined and median-filtered and run through iraf/daophot to 
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Figure 7. Spectrum of a globular cluster of NGC 4649 as seen from the three multi-object spectrographs used in this paper. Red is DEIMOS, blue is GMOS 
and green is Hectospec. Spectra are shifted to the rest frame wavelength. An arbitrary vertical offset is applied. Important absorption lines are labelled. Note 
that the Hectospec spectrum was corrected for telluric absorption, while the DEIMOS spectrum was not, which explains why the strong water absorption 
feature at ~7700A appears in the latter but not the former. The globular cluster’s ID is D28. Its apparent magnitude is 2 : = 20.3, and its radial velocity is 
1160 ± 7km s-l. 


obtain photometry for point-like objects. The instrumental photom¬ 
etry was calibrated using F555W and F814W //5rAVFPC2 pho¬ 
tometry. A final sample of globular cluster candidates was obtained 
using magnitude (9 < V < 22.5 mag) and colour (0.75 mag < 
V — 7 < 1.4 mag) cuts. Extended sources were eliminated by eye, 
leaving 250 globular cluster candidates. A GMOS multi-slit mask 
was made for Field 1 using the standard Gemini GMMPS software, 
containing 39 object slits. 

For the Fields 2-3-4-5 (program ID GN-2010A-Q-37), GMOS 
images in the SDSS g'r'i' filters were taken in 2007A and 2009A to 
select GC candidates for multi-slit spectroscopy. Colour and mag¬ 
nitude cuts and image classification using SExtractor were used to 
obtain a final sample of GC candidates in each of the four fields. 
The GMMPS software was used to create multi-slit masks, and 
with the use of filters to restrict wavelength coverage, we were able 
to put slits on between 48 to 69 globular cluster candidates in the 
four fields. Some marginal objects were included to fill the mask, 
and we also placed some slits on the diffuse galaxy light. 

Spectroscopy for Field 1 (program ID GN-2003A-Q-22) was 
obtained using GMOS-N in multi-slit mode in May and June 2003. 
A total of 8 X 1800 sec exposures were obtained with a central 
wavelength of 5000 A, and 8 x 1800 sec were obtained with a 
central wavelength of 5050 A, giving a total of 8 hours on-source 
time. We used the B600 grating, giving a dispersion of 0.45 A/pixel, 
a resolution of ~ 5.5 A, which corresponds to 320 km s~^ velocity 
resolution. Spectral coverage was typically from 3300—5900 A, 
and the seeing ranged from 0.65 to 0.9 "over the four nights of 
observation. Bias frames, dome flat-fields and CuAr arcs were also 
taken for calibration. 

Spectroscopy for the remaining four fields was obtained at 
GMOS-S in the spring of 2010. We obtained between 5 to 7.5 hours 
on-source per field, and conditions were good with seeing between 
0.5 and 1.0 arcsec. 

Data reduction was carried out using the Gemini/GMOS IRAF 
package, and consisted of the following steps: (i) a final spectro¬ 


scopic flat-field was created; (ii) object frames were bias-subtracted 
and flat-fielded; (iii) the arc frame was reduced and used to estab¬ 
lish the wavelength calibration; (iv) object frames were wavelength 
calibrated and sky-subtracted; (vi) ID spectra were extracted from 
each object frame and median combined. 

Iraf/fxcor was used to obtain radial velocities for each 
object via cross-correlation, using the MILES spectral library 
dSanchez-Blazquez et al.ll^06h as template. In this case, the tem¬ 
plate with the highest cross-correlation coefficient was used, al¬ 
though the relative scatter between templates was ~ 4 — 5 km s“^, 
well within the measurement errors. 

Before accounting for overlaps, there were ~ 185 objects with 
reliable velocities in fields 2-3- 4-5. Astrome t ry and photometry for 
this sample were obtained from Faifer et ^ (1201 ih . Globular clus¬ 
ters from Field 1 dBridges et al.l 20061) were then added to this list, 
and objects observed in more than one field were identified. 

5.3 MMT/Hectospec 

Hectospec is a multi-fibre spectrograp h mounted on the 6.5 meter 
MMT telescope dFabricant et alj|200^ . It has a 0.78 deg^ field-of- 
view which can be filled with up to 300 fibers. 

GC candidates were selected from the MegaCam and ACS im¬ 
ages. Bright objects at large radii were favoured, although this nat¬ 
urally increases the contamination from Galactic stars. Due to the 
wide Hectospec field-of-view, GC systems in galaxies surrounding 
NGC 4649 (like NGC 4621 and NGC 4638) were also observed. 
Therefore, some contamination is expected from these galaxies. 

One single Hectospec field was observed on May 16 and May 
19, 2012, using a 270 1/mm grating, with a dispersion of 1.21 
Apixel, and a spectral resolution of ~ 5 A, which corresponds to a 
velocity resolution of ~ 300 km s“^. _ _ 

Raw data were reduced as described in iMink et ^ d2007t) . 
Heliocentric radial velocities were measured by cross-correlating 
the science spectra with a Hectospec template of an M31 GC. Er- 
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rors were estimated through Monte Carlo simulations with an addi- 
tional 8 km s~^to account for the wavelength calibration (see also 
IStrader et alJ201 ih . 


6 OUTCOME OF THE SPECTROSCOPIC FOLLOW-UP 

We observed ~ 600 objects with DEIMOS, targeting ~ 500 pho¬ 
tometrically selected GCs. We found 304 spectra (including du¬ 
plicates between masks) consistent with the systemic velocity of 
NGC 4649. 23 additional GCs were discarded as “marginal” be¬ 
cause the team could not reach a consensus on whether or not the 
features seen in their spectra were actual spectral lines. In addition, 
we found 13 Galactic stars and 56 background galaxies. The latter 
were identified via the strong emission lines in their spectra. How¬ 
ever we do not measure the redshift of background galaxies. The 
remaining slits were unclassified, because the signal-to-noise was 
too low to retrieve any physical information from the spectrum. Af¬ 
ter accounting for repeated objects, the median velocity uncertainty 
of the DEIMOS GC sample is An = 11 km s“^. 

We observed a total of 187 objects with Gemini/GMOS, 
finding 172 extragalactic GCs, 14 Galactic stars and 1 back- 
gr ound galaxy. The 1 72 GCs also include the 38 GCs observed 
by IPierce et"^ l l2006l) . CMOS maps the innermost 6 arcmin of 
NGC 4649 and it overlaps nicely with the DEIMOS field-of- 
view. The median velocity uncertainty of this sample is An = 
12 km s“^. 

We observed 478 objects with Hectospec, and we found 102 
extragalactic GCs (including duplicates), 157 background galaxies 
and 119 Galactic stars. The remaining fibres were unclassified. The 
median velocity uncertainty of this sample is An = 14 km s“^, but 
in Section|7]we will explain that the effective uncertainty is larger 
than this value. 


7 REPEATED MEASUREMENTS 

We quantify how the radial velocity of a given GC compares to 
the radial velocity of the same GC observed over multiple spec¬ 
trographs. We denote radial velocities from DEIMOS, GMOS and 
Hectospec wih no, vg and nn respectively. The uncertainties of 
the repeated objects are the sum in quadrature of the uncertainties 
of the single objects. The comparison between the three catalogues 
is illustrated in Figure[8]and discussed in the followings. 


7.1 Duplicates from the same spectrograph 

We find 26 and 33 GCs repeated across different DEIMOS and 
GMOS masks, respectively. For DEIMOS, the weighted-mean 
of the velocity difference of the objects in common is (dd — 
dd) = —0.6km s“^, with an root-mean-square difference rms= 
10 km s“^. We conclude that the DEIMOS radial velocities are ro¬ 
bust, and no additional velocity offset is needed. 

For GMOS, we find field to field offsets in the range |(nG — 
iig}! = 25 — 40km s“^, which probably reflects the different in¬ 
strument set-ups of the five GMOS fields. The repeated radial ve¬ 
locity measurements from the five fields are generally consistent 
within 2a. We decided to simply average the repeated GMOS GCs 
with no rigid offset, and use the comparison with the DEIMOS 
sample to renormalize the velocity uncertainties of the GMOS sam¬ 
ple (see 117.3b . 



19 20 21 22 23 

g I mag] 


Figure 8. Repeated radial velocity measurements. Each panel shows the 
velocity difference for the objects in common between different catalogues 
(labelled on the top left) vs the magnitude in the g band. The top panel 
shows the comparison between GCs observed over different DEIMOS 
masks. The rms scatter in units of km s“^ of the velocity difference is re¬ 
ported on the top right of each panel. The velocity agreement is usually good 
(always < 100kms“^), although the velocity uncertainties from GMOS 
and Hectospec are probably underestimated (see E3)- 

7.2 Duplicates from different spectrographs 

There exist 67 GCs in common between the DEIMOS and GMOS 
catalogues. This sample has (vd — vg) = —0.6km s“^, and the 
root-mean-square is rms= 23km s“^. 

The DEIMOS and Hectospec samples have 21 GCs in com¬ 
mon. The average velocity difference is (ud — fn) = 9.3 km s“^, 
with a scatter of rms= 30 km s“^. 

Lastly, we compare the GMOS and Hectospec catalogues. 
In this case the average velocity difference is (dh — vd) ~ 
—6km s“^, with rms= 43km s“^. In the bottom panel in Figure 
[^it is clear that some of the 25 objects in common between these 
two catalogues, scatter up to 5 ct from each other. 

7.3 Calibration of velocity uncertainties 

By comparing radial velocities from different instruments we have 
seen that the average velocity difference is generally small (< 
15 km s“^), comparable to our DEIMOS uncertainties. However, 
the large scatter in (Av) compared to the small error bars of the sin¬ 
gle datapoints, implies that the velocity uncertainties of GMOS and 
Hectospec catalogues are probably underestimated with respect to 
those of DEIMOS. Therefore, we need to re-normalize the velocity 
uncertainties of the first two sp ectrographs to those of DEIMOS. 
We follow the method of IStrader et alj ( 1201 ih . and we mul- 
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tiply the uncertainties of the CMOS and Hectospec by a factor r 
until the difference between the velocities of the two catalogues 
is x^ ~ 1. This normalization counts as a first-order correction to 
the velocity inconsistencies discussed above. We find that GMOS 
uncertainties need to be inflated by a factor of r = 1.6 with re¬ 
spect to DEIMOS. Once the uncertainties have been normalized, 
we weight-average the radial velocities of the GCs in common and 
we add the unique GMOS GCs to the DEIMOS catalogue. No rigid 
velocity offset has been applied to the GMOS velocities, as this is 
negligible ((vd — vd) = —0.6 km s“^). 

Next, we compare the combined DEIMOS+GMOS catalogue 
to Hectospec. We renormalize the uncertainties of the Hectospec 
GCs by a factor of r = 2.4. Also in this case the velocity difference 
between the two datasets is negligible Av = 5 km s“^, therefore 
no offset is applied. We merge Hectospec to the master catalogue 
following the above procedure. Our master catalogue is composed 
of 447 spectroscopically confirmed stellar objects (GCs, UCDs and 
dwarf galaxies). 


8 BUILDING UP THE SPECTROSCOPIC MASTER 
CATALOGUE 

Before proceeding, we assign photometric measurements (i.e. mag¬ 
nitudes and sizes) to all confirmed objects. We match the position of 
confirmed GCs with all sources, both point-source and extended, 
extracted from ground based and HST imaging. We transform the 
ground based colours of the confirmed GCs from {g — i) to {g — z) 
using Equation|2 This is more convenient than transforming (g—z) 
to (g — i) because 60 per cent of the confirmed GCs have ACS pho¬ 
tometry. 


8.1 Clipping GCs from neighbour galaxies 

The large field-of-view of Hectospec allows us to confirm GCs 
in some galaxies near NGC 4649 (see Figure [TJ. These GCs 
are treated as interlopers because they do not follow the grav¬ 
itational potential of NGC 46 49. The galaxies whi ch may con¬ 
taminate our sample are (from l&ppellari et al.ll201 ih : NGC 4621 
(usys = 467km s-1), NGC 4638 (usya = 1152km s’^), 
NGC 4647 (Uaya = 1409 km s-\ from NED), NGC 4660 (Usys = 
1087km s“^). We remind the reader that NGC 4649 has Usys = 
1110kms“^ dCattpellari et 

Of the 7 GCs around NGC 4621, five are probably associated 
with this galaxy because their radial velocities are very similar to 
that of this galaxy as we will show in Section [9l The remaining 
two GCs (those in the NGC 4649 direction) have radial velocities 
consistent with NGC 4649 and they probably belong to this galaxy. 
Using geometrical criteria, we find 4 GCs associated with the ellip¬ 
tical galaxy NGC 4660. 

We flag one GC as a member of NGC 4638 because it is 
very close to the center of this galaxy. However, the membership 
of this GC is impossible to confirm because the systemic velocity 
of NGC 4638 is almost identical to that NGC 4649. 

GCs belonging to the spiral galaxy NGC 4647 (2.6 arcmin 
from NGC 4649) m ay also contaminate our GC sample. However, 
I Strader et alj j2012h estimated that only ~ 10 objects in their GC 
catalogue are actually associated with NGC 4647. From our spec¬ 
troscopic GC catalogue, we exclude one GC because of its prox¬ 
imity to NGC 4647 and because the spiral arms compromise the 
photometry of this object. 


8.2 Clipping surrounding galaxies 

In addition to GCs, we also measured radial velocities for some 
dwarf galaxies in the field. 

With DEIMOS and GMOS, we observed the dE VCC 1982 
[v = 945 ± 4km s“^) and the nucleated dE VCC 1963 (u = 
1090 ± 13 km s“^). For comparison. lUee et'H] ( l2008h found ve¬ 
locities of 976 ± 42km s“^ and 1027 ± 53 km s“^, respectively. 
We also obtained a radial velocity for a very lo w surface brightness 
nucleated dEs (ID 81 from ISabatini et al.l2005h . We found its radial 
velocity to be u = 1186 ± 22 km s“^. 

With Hectospec, we obs erved two dwarf g alaxies present in 
the photometric catalogue of iGrant et al.l feOOsI) . The first, is the 
nucleated dE VCC 1951, for which we findn = 1089±32 km s“^, 
consistent with the measurement from SDSS DRIO v = 1106 ± 
13 km s“^. The second, is the dE VCC 1986, for which we measure 
the redshift for the first time, finding v = 850 ± 37km s“^. 

The object M60-U CD1 is also treated as a galaxy and ex¬ 
cluded from our sample jStrader et al.ll2013h . 

One GC is only 7 arcsec away from a faint dwarf galaxy with¬ 
out spectroscopic confirmation, but we assume this is a projection 
effect and that these two objects are not physically connected. The 
same assumption was adopted for another GC which is ~ 8 arcsec 
away from the centre of a small galaxy without a redshift measure¬ 
ment (J124334.56-f 112727.3 from SDSS DR8). 


8.3 Additional GC datasets 


We searched the SDSS archive for point sources with spectro¬ 
scopic confirmation within 1 deg from NGC 4649. We found 
one unique object, which we added to our master catalogue, 
consistent with being a GC associated with NGC 4649 (ID 
SDSSJ124349.56-f 113810.3). We found four other GCs with 
SDSS spectroscopic confirmation, but these were not added to our 
catalogue because they already observed and confirmed with Hec- 
tospeC; _ 

lUee et al] l l2008l) published radial velocities for 93 GCs us¬ 
ing CFHT/MOS with a median velocity uncertainty of 50 km s“^. 
Their observations explored an area « 14 x 14 arcmin^, which 
overlaps the spectroscopic observations from our paper. In fact, we 
found 61 GCs in common between our spectroscopic master cata¬ 
logue and Lee et al.’s dataset. The median velocity difference is 38 
km s“^, with an rms of rms=156 km s“^. A more careful inspec¬ 
tion of the dataset, shown in Figure [O] reveals that the large rms 
value is driven by two objects (ID 201 and ID 226 using Lee et 
al. nomenclature) which scatter up to 600 km s~^ or (Sa) from the 
ideal An = 0km s“^ line. Removing these two extreme objects, 
the velocity difference rms decreases to 113 km s“^, still a factor 
of two larger than measured in our dataset. 

The magnitude o f this effect is compa rable to that of the out¬ 
liers found in M 87 by Strader_e£^ j201lj) when comparing to the 
CFHT/MOS dataset of Hanes et ^ 200 ih . We investigated the rea¬ 


son for this offset. We excluded mismatch because the outliers are 
not in overcrowded regions. We excluded a low signal-to-noise is¬ 
sue, because the two objects have intermediate magnitudes and be¬ 
cause fainter duplicates have radial velocities in reasonable agree¬ 
ment with our dataset. As we could not identify the reason for this 
disagreement, we decided not to include the Lee et al. dataset in our 
DEIMOS catalogue because additional unknown outliers amongst 
the 32 unique Lee et al. GCs might skew our GC kinematic results. 
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Figure 9. Comparison with the data of iLee et"^] i20Q8h . The plot is the 
same as Figure^ but for our com bined spectroscop ic master catalogue ver¬ 
sus the CFHT/MOS GC dataset of iLee et'H] i200^ . Two objects scatter up 
to 700km s~^ with respect to the Av = 0 dashed line. 
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9 OVERVIEW OF THE SPECTROSCOPIC CATALOGUE 

We compile a list of 448 unique spectroscopically confirmed ob¬ 
jects around NGC 4649. After clipping dwarf galaxies and GCs 
from neighbour galaxies, we are left with 431 unique GCs belong¬ 
ing to NGC 4649. The spectroscopic catalogue is given in Table 
[2l In the following, we review different projections of our param¬ 
eter space, which involves photometric quantities (magnitude and 
colours) as well as physical quantities, such as the radial velocity 
and the physical size. 

9.1 Spatial distribution 

We start by looking at the positions on the sky of all confirmed ob¬ 
jects (Figure [Tot. We confirm GCs up to 100 kpc from NGC 4649. 
DEIMOS and GMOS contribute to the GCs in the innermost 30 
kpc, whereas all objects outside this radius are from Hectospec. In 
the outermost regions, mostly blue GCs are confirmed as expected 
from the shape of the GC surface density profile. 

From Figure [To] is also clear that the GC systems of nearby 
galaxies are spatially separated from NGC 4649. Additional con¬ 
tamination from these galaxies is unlikely. The paucity of data 
in the immediate North-West of NGC 4649 is due to NGC 4647, 
which prevents us from detecting any source. 


9.2 Globular cluster velocity distribution 

In Figure [TT] we plot the GC velocity distribution vs. galactocen- 
tric distance. This diagram is useful for identifying outliers which 
may scatter from the GC system velocity distribution. The latter is 
shown in the right panel of Figure [TT] The skewness and kurtosis 
of the GC velocity distribution (excluding GCs from other galax¬ 
ies) are s = —0.05 ±0.11 and ac — 3 = —0.21 ±0.23, respectively, 
suggesting Gaussianity. This is reassuring because we will invoke 
Gaussianity later on. If we consider only GCs within the innermost 
100 arcsec, we estimate the systemic velocity of the GC system to 
be Vsys = 1115 ± 28 km s~^, very consistent with literature values 
for NGC 4649. 

We bin the GCs by radius and calculate a, which is the stan¬ 
dard deviation of the binned velocities with respect to the systemic 
velocity of the galaxy. This allows us to draw a 3a envelope as a 
tool to flag GCs deviating more than 3a from the local GC velocity 
distribution. As can be seen in FigurefTT! none of the objects shows 
such a deviation.. 

The velocity distribution of GCs is distinct from that of Galac¬ 
tic stars and from background galaxies. Therefore, a velocity cut at 



Figure 10. Spatial distribution of the spectroscopically confirmed objects. 
NGC 4649 is at (0, 0). Top and right axes show the co-moving physical 
scale at the distance of NGC 4649. Confirmed GCs are divided into blue 
and red subpopulations and coloured accordingly. The green star marks the 
position of NGC 4647. Orange boxes are dwarf galaxies discussed in the 
text. Ellipses are schematic isophotes at 3,6, 9 Re- Position of some of the 
surrounding galaxies and their confirmed GCs are marked with coloured 
stars and squares, respectively. North is up, east is left. 


V = 400 km ensures no contamination from Galactic stars. We 
note that for R > 700arcsec, the majority of GCs have blue-shifted 
radial velocities. Within R = 200 arcsec, neither Galactic stars nor 
background galaxies were found thanks to the superior spatial res¬ 
olution of ACS.. Confirmed GCs from other galaxies are clustered 
around the systemic velocity of the host galaxy. 


9.3 Trends witb colour 

Figure [T2 illustrates how some physical and photometric proper¬ 
ties correlate with the observed GC colour. The quantities that we 
explore are the magnitude, the physical size and the radial ve¬ 
locity V. 

In the top panel we plot the colour histogram of all the HST 
candidates compared to the outcome of our spectroscopic follow¬ 
up. The overall colour distribution is well sampled with confirmed 
GCs. Bimodality is preserved, albeit there may be a deficit of redder 
GCs. 

The correlation of colours with velocity (Panel b) suggests that 
the blue GCs scatter with respect to Pays more than the red GCs. In 
other words, the velocity dispersion is larger than that of the red 
GCs. We will quantify this statement in the next section. 

Panel c shows how colour correlates with size. The average 
GC size of the sample is rj, « 2 pc. We confirm 10 objects with 
Th > 10 pc. We note that the catalogue of ! Strader et al.l l l2012h lists 
22 objects with rn > 10 pc, therefore our return rate is about 50 
per cent. We confirm all three o bjects with > 20 pc, two of 
which were already confirmed in lLee et^ ( l2008l) . The paucity of 
UCDs in this galaxy prevents us from any kinematic analysis of the 
UCD sample. We define UCDs as objects brighter than uj Cen and 
with Vh > 10 pc, w h ere we have adopted the size cut adopted by 
iBrodie eT ^ ( I2OIII) or lZhang et alj ( l2()15h . 
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ID 

ID(S-|-12) 

RA 

DEC 

{9 - 2) 

Z 

R 

V 

5v 



Ideg] 

Ideg] 

[mag] 

[mag] 

[arcsec] 

[km/s] 

[km/s] 

412 

- 

191.09555 

11.45407 

1.095 

20.177 

694 

1240 

32 

295 

- 

190.85417 

11.52163 

0.991 

20.455 

237 

1034 

23 

114 

C4 

190.94971 

11.59446 

1.140 

19.943 

197 

1247 

22 

59 

B12 

190.92905 

11.53665 

1.574 

19.759 

75 

612 

6 

339 

- 

191.04245 

11.57867 

0.948 

18.483 

420 

1450 

28 

353 

- 

191.01315 

11.61024 

1.331 

21.138 

387 

771 

14 

96 

B87 

190.97552 

11.51810 

1.361 

21.552 

233 

1076 

21 

160 

D28 

190.90526 

11.52134 

0.948 

20.352 

129 

1160 

7 

333 

- 

190.80938 

11.57225 

1.529 

21.760 

354 

1424 

8 

295 

■ 

190.85417 

11.52163 

0.991 

20.455 

237 

1034 

23 


Table 2. Catalogue of the 431 spectroscopically confirmed GCs associated with NCG 4649. ID(S+12) is the ID of the GCs also present in Istrader et alJ 
Right ascension and declination are given in degrees, {g — z) and 2 : are either from Strader et al., or have been converted from gri magnitudes as 
described in El R is the galactocentric radius, v is the observed radial velocity of the objects. The quoted uncertainties on v, 5v, have been re-normalized to 
the DEIMOS uncertainties (see text). The full table is available online. 
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Figure 11. GC velocity distribution vs. galactocentric radius. Shown are all 447 spectroscopically confirmed objects. Left panel: observed radial velocities vs 
galactocentric distance R, shown in physical units on the top axis. Error bars are the velocity measurement errors re-normalized to the DEIMOS uncertainties 
(see text for details). Grey points ai'e objects associated with NGC 4649; orange squai'es are the dwarf galaxies discussed in the open black circle at 
R ~ 100 arcsec is the object M60-UCD1. Cyan points are GCs associated with NGC 4621; magenta points are GCs associated with NGC 4660; the brown 
point is a GC probably associated to NGC 4638. Open boxes show the position in this diagram of the three galaxies above. Green stars are Galactic stars. The 
±3 (t envelope of the NGC 4649 GCs (grey points only) is shown as a blue lines. Horizontal lines show v = 0 km s“^ and v = Vsys = 1110 km s“^, which 
is the adopted systemic velocity of NGC 4649. Vertical thin black lines show the position of confirmed background galaxies with arbitrary velocities. Right 
panel: velocity distribution histogram for GCs/UCDs belonging to NGC 4649 only (grey histogram), and Galactic stars (green histogram). Both histograms 
have the same bin-step of 100km s“^. GCsAJCDs are well separated from Galactic stars and background galaxies at all radii. 


It is interesting that very extended objects are either the red¬ 
dest, or the bluest. Among the blue objects, only 2 can be defined 
as “classic” UCDs, because along with being larger than GCs, they 
are also bright (Mz < —11.0 mag). The remaining blue objects 
ha ve intermediate ma gnitudes and are the faint UCDs discussed 
in lSrodie et al.l ( 1201 ll) . Regarding the confirmed extended objects 
with red colours (four in total), these are generally fainter than 
z > 21, and never larger tha n rn ~ 20 pc. These ob jects are more 
likely linked to faint-fuzzies dBrodie & Larsen|[2002h . 

Lastly, the bottom panel shows the colour-magnitude diagram 


of both candidate and confirmed GCs. We notice that some are 
brighter than cj Cen (see also Figure O- The M60-UCD1 object 
has z: = 15.6 mag and is not visible in the plot. The blue-tilt is 
clearly visible, as well as the fact that the red subpopulation has a 
larger colour spread compared to the blue subpopulation. We note 
the presence of ~ 10 objects with z < 20 mag and unusually blue 
colours {g — z) < 0.95 mag, which are found outside the ACS 
field-of-view and therefore lack of size information. It will be in¬ 
teresting to determine why these objects do not tilt towards redder 
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Figure 12. GC trends with colour. Each panel shows the {g — z) colours 
of the spectroscopically confirmed GCs (large grey points) and of the HST 
GC catalogue (small black points). Coloured points are confirmed GCs from 
other galaxies as defined in Figure lTTI The vertical line marks the sepai'ation 
between blue (on the left) and red GCs (on the right). Objects with ground- 
based photometry lack size measurements, and they were given = 100 
pc, as can be seen in Panel (c). The black and grey histograms in Panel (a) 
are the HST dataset and the confinned GCs, respectively; the horizontal line 
in Panel (c) separates UCDs from GCs; the horizontal lines in Panel (d) ai'e 
the magnitude of uj Cen at the distance of NGC 4649 and the turn-over- 
magnitude. 


ISchuberth et al.|[2010al : IPota et al.ll2Q13h . We now quantify these 
findings and we derive physical quantities capable of characteriz¬ 
ing the kinematics of the NGC 4649 GC system. We are interested 
in its rotation amplitude Urot, velocity dispersion a, kinematic po¬ 
sition angle 6o, and root-mean-square velocity dispersion Urms- 


10.1 Method 


We start by calculating the root-mean-square velocity dispersion 
profile Urms, which measures the standard deviation of GC radial 
velocities from the systemic velocity of the galaxy: 


V 


2 

rms 



z 




(3) 


where Awi is the measurement uncertainty of the i-t h GC. Uncer¬ 
taintie s on frms are obtained using the formula of IPanese et alJ 
( Il980h . The Urms can be thought of a measure of the specific ki¬ 
netic energy of a system Urms ~ frot + and does not require 
any assumption about the detailed velocity distribution of the sam¬ 
ple. However, we are also interested in the relative weights of Wrot 
and o in the final value of Urms, which cannot be deduced from 
equation]^ 

Therefore, in order to estimate a, Vrot and PAkin, we adopt 
a model for the GC kinematics. This approach comes at the cost 
of assuming that the line-of-sight velocity distribution is Gaussian, 
but we have shown that this is a reasonable assumption to a first 
approximation. __ 

W e use the rotation model dCote et al.l 1200 ll : IProctor et al.l 
l2009li : 


ttmod — Vsys i 


ttrot 



^ tan(Pyt-PAkin) ^^ 


(4) 


where Wrot is maximum amplitude of the sine curve; a is the stan¬ 
dard deviation of the data points with respect to Vmod; the signs 
H- and — are applied if the PA lies inside or outside the range 
{9 — 9o) = [—7r/2, -|-7r/2], respectively; q = 0.84 is the axial ratio 
of the GC system, assumed to be the same as that of NGC 4649. 

To find the set of (urot, cr, PAkin) that best reproduces the 
data, we minimize the log-likelihood fun ction (i.e. the function) 
teergond et aDl2006l ; [^der et alj20llh : 


X 


2 




-b ln[o-^ -b (Aui)^] 


(5) 


where a is the standard deviation of the datapoints with respect to 
the model. 

Uncertainties are derived with “bootstrapping”. We randomize 
the sample of interest 1000 times and we fit the data with equations 
|4]and|3each time. We then infer the 68 per cent (la) error from 
the cumulative distribution function of each free parameter. 


colours like to bulk of blue GCs. We defer such analysis to a future 
paper. 


10 KINEMATICS 

So far, we have qualitatively discussed some properties of the 
spectroscopic sample, with no assumptions about the shape of the 
GC velocity distribution. This has suggested that blue and red 
GCs have different properties, as found in other GC systems (e.g.. 


10.2 Independent tracers: stars and planetary nebulae 

It is interesting to compare our GC kinematics results with those 
from two independent kinematic tracers: diffuse s tarlight from Fos¬ 
ter et al. (2015) and planetary nebulae (PNe, iTeodorescu et al.l 

l201lh . 

The stellar data extend out to 2Re. We do not compare our 
results with long-slit or integral-field-unit data because these extend 
out to ~ 0.5i?e and they do not overlap with our GC data in the 
innermost regions. 
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The PNe dataset consists of 298 objects. The PN that is fur¬ 
thest from the galaxy centre is at 400 arcsec, which corresponds to 
the radial extent of the DEIMOS GC sample. We infer the kine- 
ma tics of PNe using equatio ns [3] and as done for GCs. Follow¬ 
ing |^odorescu|et^ ( l201lh . an average measurement uncertainty 
Au = 20 km s“^ is assumed for all PNe. 

The kinematics and sp atial distributions of s t ars and PNe have 
been shown to be similar dCoccato et al.l l2009l : iNapolitano et al.l 
I 2 OO 9 I 1 because they are different manifestations of the same stel¬ 
lar population. It is less clear how the properties of PN systems 
general ly compare with those of the GC system within the same 


galaxy 1 WpodleveLal ][20Ki yFoster et ai]|201 iLlporbes et alj2012l : 

IPota et al.ll2013l:ICoccato et al.ll2013lh 


11 KINEMATIC RESULTS 

We divide our sample into two subsamples. The “GC” sample is 
made by objects with i > 20 mag. We also require sizes < 10 
pc for the blue GCs ((p — z) < 1.2), if available, and ru < 20 
pc for the red GCs in order to include faint-fuzzies as part of the 
red GC sample. This choice is motivated by the findings that faint 
fuzzies are associa ted with galaxy disks dChies-Santos et ai]|2013l : 
iForbes et al.ll20l4l . although it is unclear if NGC 4649 possesses 
an embedded disk, as we discuss in Section If physical sizes 
are not available, the object is still included in the GC sample if 
the magnitude criterion is satisfied. This GC sample includes clas¬ 
sic GCs fainter than u Cen, reflecting our attempt to prune from 
the sample extended objects without size measurements, but it may 
also include faint extended objects with no size measurements. The 
“bright” sample includes objects with 2 < 20 mag with no con¬ 
straints on the size, and represents bright GCs and UCDs. 

Fach group is divided into blue and red GCs adopting a con¬ 
servative colour cut at {g — z) = 1.2 mag. This choice introduces 
mixing between the two GC subpopulations for (g — 2 ) ~ 1.2±0.1 
mag. Therefore, we also investigate the effects of the colour cut on 
the results by selecting blue GCs with {g — z) < 1.1 mag and red 
GCs with (g — z) > 1.3 mag, respectively. 

We start with a sample of 431 GCs. We exclude one isolated 
GC at 1500 arcsec from NGC 4649, and one GC because of its 
proximity to a bright star, which affects the magnitude measure¬ 
ments. In summary, we analyse the kinematics of a sample of 429 
objects around NGC 4649. 


11.1 The rms dispersion profile 

We compute the GC Urms velocity dispersion profile in radial bins 
wifh roughly 50 objects per bin. For the stars, we compute u^ms = 
Vrot + using the tabulated values from Foster et al. (2015). 

We calculate Urms for GCs and PNe in rolling radial bins with 
50 objects per bin. The results are shown in Figure [T^ The Urms 
profiles from the three tracers are diverse. The red GC Vrms pro¬ 
file decreases monotonically with radius, whereas the Vrras profile 
of blue GCs and PNe are bumpy and akin to each other, which is 
unexpected when compared to other galaxies. 

We also note that the PNe results agree with the stellar Vrms 
profile in the region of overlap. On the other hand, the results for the 
red GCs seem to be systematically larger compared to the stars. It is 
unclear if this effect is real, or is du e to unknown syste matics in the 
spectroscopic measurements (e.g.. lArnold et al.ll201lh . This result 
is still significant when a more extreme GC colour cut is adopted. 



R [arcsec] 

Figure 13. Root-mean-square velocity profile of kinematic tracers in 
NGC 4649. Shown as think lines are the rolling-binned Urms profiles (equa- 
tion[3 for blue and red GCs coloured accordingly, respectively. PNe results 
are shown in green. Stellar data are shown in purple. Dashed lines are Icr 
envelopes. 


11.2 Results from the kinematic modelling 

The best fit parameters to eq. |4]for GCs and PNe sample (without 
any radial binning) are given in Table GCs from surrounding 
galaxies are discarded from the analysis. 


11.2.1 Kinematic radial profiles 

Radial kinematic profiles for GCs, PNe and slars, are shown in Fig¬ 
ure [14] We show running average profiles, wifh ~ 40 GCs per bin. 
Focusing on fhe red GCs (lefl panels), this sample shows a nearly 
constant rotation with radius. The rotation axis is very consistent 
with the photometric major axis of the galaxy (PAkin ~ 270 deg). 
The velocity dispersion is bumpy and decreases with radius. The 
dip in rotation velocity between 100 and 200 arcsec is intriguing 
and it will require future in-depth analysis. One possibility is that 
this feature is spurious because it occurs in the radial range where 
the spiral galaxy NGC 4647 prevents us from detecting GCs, there¬ 
fore biasing the kinematic fit (see also Figure do). Alternatively, 
the velocity dip may be real and linked to the recent inte raction be¬ 
tween NGC 4649 and the galaxy remnant M60-UCD1 dSeth et al.l 
l2014h . which lies at about 90 arcsec (7.2 kpc) along the minor axis 
ofNGC 4649. 

The blue GCs (middle panels) also show hints of rotation at all 
radii. The rotation amplitude is lower with respect to the red GCs. 
As with the red GCs, the rotation direction is generally consistent 
with the major axis of the galaxy. Minor axis counter-rotation is 
also detected at 300 arcsec. The velocity dispersion is constant at 
(T = 240 ± 30km s“^, and higher than that of the red GCs, as 
found in other galaxies. The GC kinematic results are generally in¬ 
dependent of the colour cut, as shown by the open points in Figure 

□3 

Lastly, the PNe (right panels) show a nearly constant rotation 
pr ofile, but the kinematic position angle twists with radius, as found 
bv ICoccato et al.l ( 1201 3tl . The velocity dispersion profile overlaps 
nicely with the stellar results in the innermost regions. We see a 
cr-bump at 100 arcsec, and a declining velocity dispersion profile 
intermediate between the results for the blue and red GCs. 

Figure [T3 shows that the three kinematic tracers are similar 
to each other in some aspects, e.g., rotation axis and rotation am¬ 
plitude, but their velocity dispersions behave very differently. The 
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Figure 14. Kinematic radial profiles. Red, blue and green lines represent red GCs, blue GCs and PNe respectively. Open circles are the results for the reddest 
and the bluest GCs, respectively (see text). Grey points are stellar data from Foster et al. (2015). Each panel represents a solution to equation [4] for a given 
kinematic tracers. From the top to the bottom we show the rotation velocity Vrot, the velocity dispersion cr with respect to Vrot, and the kinematic position 
angle PAkin. Rotation occurs near the major axis at PA ~ 270 deg (marked with a black line). The dashed line indicates rotation along the minor axis. 
Running averages are performed with ~ 35, 40, 45 objects bin for red GCs, PNe and blue GCs, respectively. The co-moving step is 5 objects. Rotation along 
the major axis is detected at all radii. 


sample size of the bright group is too low to compute radial pro¬ 
files. However, the azimuthally averaged quantities from Table [3] 
indicate no significant difference between this group and the results 
discussed above. 


tation is suggested for a group of blue GCs with (g — z) ~ 0.9±0.1 
mag, as detected at intermediate radii (Figure [T4t . 


11.2.2 Kinematic colour profiles 

We have shown that the kinematic modelling results are not very 
sensitive to the assumed colour cut. We now bin the data by colour, 
rather than by galactocentric radius, and measure the azimuthally 
averaged kinematics of the objects within a certain colour bin. This 
exercise is independent of the colour cut. The results are shown in 
Figure [TS] in the form of co-moving bins with 50 objects per bin. 
Comparison with the colour GC gradients in Figure [^ and Figure 
I14lis informative. 

The results in Figure ED show a double-peaked rotation pat¬ 
tern corresponding to the rotation of the blue and red GCs. The 
central panel shows a declining velocity dispersion in the range 
0.8 < {g—z) <1.0 mag, which flattens at redder colours. More in¬ 
terestingly, the GCs with {g — z) > 1.55 mag have a large velocity 
dispersion and negligible rotation, because they are very centrally 
concentrated (see Figure[^. The rotation is consistent with the pho¬ 
tometric major axis of the galaxy. On the other hand, minor axis ro- 


12 DISCUSSION 


NGC 4649 is a normal looking elliptical galaxy with a relatively 
normal GC system. From a photometric point of view, we have 
shown that the GC system of NGC 4649 presents typical features 
seen in other elliptical galaxies of similar mass. We confirm statisti¬ 
cally signifi cant GC colour bimodality, as found in earlier work on 
NGC 4649 jBridges et al.ll2006l:IPierce et alj2006l:iree et alj2008l: 


iFaifer et^ 201 ll : Chies-Santos et al. 201 ll). The G C spec ific fre¬ 

quency within 50 kpc Sn = 3.5 ± 0.1 l lFaifer et aljEoilh is con¬ 
sisten t with values obtained for galaxies of similar mass ( Peng et aP 
I 2 OO 8 I) . The two GC subpopulations have different spatial distribu¬ 
tions (with the blue GCs being more extended than the red GCs) 
and negative colour gradients in the innermost regions. 

From a kinematic point of view, the GC system of NGC 4649 
is less normal when compared to other systems, because of the sig¬ 
nificant rotation detected at all radii for both GC subpopulations, 
which is a rare occurrence in elliptical galaxies of this mass. In what 
follows, we discuss possible formation scenarios for NGC 4649 
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Sample 

N 

7/rot 

a 

HAkin 



[km s“^] 

[km s“^] 

[degree] 


GC sample with {g 

— i) < 1.2 and [g 

•-i)> 1.2 

All 

366 


229tl 

256tl? 

Blue 

212 

Sltlo 

237t?o 

244tS 

Red 

154 

lists 

216t?o 

269tl3 


GC sample with {g 

— f) < 1.1 and {g 

•-i)> 1.3 

All 

279 

lists 

229t8 

246tS 

Blue 

170 

qq+28 

°®-28 

238tii 

23ltS 

Red 

109 

isotS 

209tii 

260tlt 

Bright sample with [g — i) < 1.2 and {g — i) > 1-2 

All 

60 

68tS 

273tS 

971 +^® 

Blue 

35 

6ltS 

26lt43 

-*^^^-68 

Red 

25 

9^0 + 55 

900 + 32 
^dO_3g 

293tl7 

PNe sample 

PNe 

298 

98tS 

234tS 

267tS 

Table 3. Results of the kinematic modelling. Listed are the best ht parame- 


ters to eq.|4]for different sub-samples, each containing N objects. 



0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

( g - z) [ mag] 


Figure 15. Kinematics as a function of GC colour. Each rolling bin contains 
50 GCs. 1(T envelopes are shown as thin lines. The solid and dotted lines 
indicate rotation along the major and minor axis, respectively. A colour cut 
at {g — z) = 1.2 mag is used to derive the kinematic radial prohles. 

based on the results from this paper and additional data from the 
literature. 

12.1 Comparison with simulations 

We start by looking at the inne rmost stellar kin ematics in 
NGC 4649. We use the tools given bv iNaab et al.l l l2014h to link the 
observed innermost stellar kinematics from ATLAS®° with simu¬ 
lated galaxies with known formation histories at z = 0. NGC 4649 


is classified as fast-rotator bv lEmsellem et al] | |201 ih . with param¬ 
eter Ae = 0.127 within IRe and ellipticity te = 0.156. New ex¬ 
tended stellar kinematics data from Foster et al. (2015), have shown 
that the value of A rises up to 0.6 at 2.5 Re, meaning that the galaxy 
becomes more and more rotation dominated in the outer regions. 
Moreover, we find that: 

• The pixel-by-pixel values of Vja (where V is the observed ra¬ 
dial velocity) is weakly anti-correlated and correlated hs and / 14 , 
respectively. The Pearson coefficients for the ATLAS®® data are 
Pv/a,h 3 = —0.08 and pv/<T,hi = 0.24. This is shown in Figure 
[T6l where we show both ATLAS®® data and DEIMOS stellar data 
from [poster et al.l ( 1201 5h within \Re. The parameters /13 and /14 
describe the shape of the line-of-sight velocity distribution within a 
given spaxel. The correlation with V/cr is indicative of axisymme- 
try or of the presence of a disk. 

• The stellar velocity dispersion is peaked at the centre, and de¬ 
creases with radius (see Figure [T4t. 

• Isophotes are ellipti cal (not disky) and there is no evidence o f an 
embedded stellar disc dKormendv et alj200^ : IVika et alj|2012h . 

According to the recipe of iNaab et alj ( l2014h . galaxies with 
the above properties formed via one late (2 < 1 ) gas-poor ma¬ 
jor merger and many more minor mergers. Major mergers here are 
meant as those with up to 1:4 mass-to-mass ratio. Less than 18 per 
cent of the stellar mass formed in-situ, owing to the relatively low 
amount of gas involved in the merger. The remaining mass frac¬ 
tion, the one which makes up most of the stellar halo, was accreted 
via minor mergers. The strength of the correlation between ( 13 , (14 
and Via can be linke d to the amount of gas involved in the merger 
dHoffman et al.ll200^ . However, the observed correlations in Fig- 
ure[T^are not strong, and therefore we cannot constrain the amount 
of gas, if any, involved in the major merger. 

In principle, the observational properties listed above are also 
compatible with a gas-rich major merger remnant, but these rem¬ 
nants are usually less massive than logM* = II.SMq and they 
have disky isophotes, w hereas NGC 4649 has logM* = II.TMq 
I Cap2ellariet_alj l2013h and no strong evidence for a stellar disk 
I Vika et alj2012h . Spectral energy distribution fitting ofNGC 4649 
implies an average stellar age of 9.5 Gyr, no rec ent star formation 
and negligible gas content dAmblard et al.l2oT^ . 

Further evidence supporting a dry major merger formation 
pathway for NGC 4649 comes from the deficit of stellar mass 
(Mdef) at the centre dLauer et alj2005l : lKormendv et al.l2009h . The 
interaction between two supermassive black holes (SMBHs) dur¬ 
ing a d ry merger can evacuate stars from the central regions of the 
galaxy dMerritll2006h . If the mass of the SMBH of the merger rem¬ 
nant is known Msmbh, the value of Mdef can be connected to the 
number N of dry mergers that a galaxy has undergone, such that 
Mdef /Msmbh ~ 0.5 per merger. IduIIo & Grahamld 20 l^ inferred 
Mdef/M smbh = 1-22 for NGC 4649, meaning that this galaxy 
underwent the equivalent of « 2.5 dry major mergers during its 
assembly, in general agreement with t he predicted merger rate of 
galaxies with comparable stellar mass dConselice & Arnolal2009l : 
lKhochfaretai]|201ll: iNaab et al.ll2014l) . 

In summary, formation via dry major mergers and accretion 
can qualitatively reconcile the central stellar kinematic profiles with 
predictions from galaxy formation simulations. We now consider 


^ Ae is a proxy for the galaxy angular momentum wit hin 1 Re and it 
is equ ivalent to the luminosity weighted ratio (nrot/cr) dEmsellem et alj 
1200 -^. 
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Figure 17. Ratio of the rotation amplitude over velocity dispersion as a function of galactocentric distance. GC data are for the classic GC sample only (i.e. 
from ^11.2.1l and Figure [HI . The left and right panels show the red GCs plus PNe (in green), and blue GCs respectively. Stell ar data ai'e shown in bo th panels 
as small grey point s. The thin dashed line and the thin solid line are the [vrot/cr) predicted in a 1:1 dissipationless merger IPi Matteo et alj|2009al) and 1:3 
dissipative merger ICretton et alj200lll . respectively. Thick dashed lines are the maximum (urot /c) predicted in the 1:1 dissipationless merger simulations of 
[Foster et ^ ^201 ih (model 4 in their Table Cl). 
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Figure 16. Pixel-by-pixel comelations between IZ/cr and /13 and . Gray 
points are ATLAS®^ data, whereas black points are DEIMOS stellar data 
within 1 Rs- 


whether this picture is also consistent with the GC constraints dis¬ 
cussed in this paper. 

As discussed above, rotation in elliptical galaxy outskirts can 
be interpreted as the dynamic al imprint of a major merger between 


two rotating di s k prog e nitors dCote et al.l200l]: 

IVitvitska et alj 120021 : ISekki et ^ 20ok 
or bet ween two pressure-supported spheroids 1 Di Matteo et al.l 
l2009ah . The rotation of GCs (and other trace particles) in galaxy 
outskirts is also expected to have increased af ter a major merger 
teekki et al.ll2005l : [McNeil-Movlan et al.ll2012h . This effect can be 
quantified through the ratio (frot/o"), which we show in Figure [77] 
for stars, GCs and PNe. 

The left panel of Figure[T7]shows that stars have an increasing 
(wrot /o) profile, whereas red GCs and PNe follow a bumpy pattern. 
We note that (urot/f) is always < 1, and that the profiles seem to 
decline outside 20 kpc, which may indicate the turning point of the 
galaxy rotation curve. On the other hand, the blue GCs (right panel) 
have less rotation compared to other tracers, with a rising (urot/cr) 
profile at large radii. 

Results for the stars are compar ed with the predicted (urot/o") 
profile of a gas-rich m ajor merger J&e tton et a lj|200lh and of a 
gas-poor major-merger foi Matteo et alj|2009ah . Both simulations 
can reproduce the increasing stellar Urot /o profile in the region of 
overlap. The predicted Vrot /o' keeps rising with radius, whereas our 
data show a more irregular rise towards higher Urot/o". 

Similarly, the Urot/o" profile of GCs is co mpared with t he dis- 
sipationless major mergers simulations from [Foster et"^ ( 1201 ih 
(thick lines in Figure [TtIi. These simulations can qualitatively re¬ 
produce the maximum Urot/o" of blue and red GCs results for a 
wide range of orbital configurations (see Table Cl in [Foster et al.l 
l201lh . Therefore, the radial increase of Vrot/o for the GC is 
qualitatively consistent with a dry major merger picture. We note 
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that kinematic predictions for GC kinematics in gas-rich ma- 
jor merger remnants are currently not available in the literature 
(lKrui |ssen et alJ2012h . 


Schauer et al. 1 12011 showed that major merger simulations 


predict that the velocity dispersion profile of the kinematic tracers 
of the remnant should be “bumpy” at intermediate radii. This is a 
permanent feature within the galaxy , and it has been observed in 
some systems te.g.. lPota et al.l 2013ll . We find that the velocity dis¬ 
persion of red GCs and PNe in NGC 4649 is bumpy, especially at 
~ 100 arcsec (see Figure m. By extrapolating the stellar veloc¬ 
ity dispersion profile at large radii, we estimated that the bumps 
have an ampli tude of ~ 50 km s~ ^, in general agreement with the 
predictions of ISchauer et al.l l l2014h . 

It is known that a galaxy’s metallicity profile can be contin¬ 
uously reshaped during its life cycle. Dissipative collaps e tends to 
create steep metallicity radial profiles dPinino et al.ll201ol) . whereas 
major mergers can steepen or flatten metallicity gradients, depend - 
ing on the initial conditions of the progenitors dpont et alj |201 lb . 
Similar r esults can also be achiev ed through dissipationless major 
mergers dOi Matteo et al.ll^09bll . Other proces ses, such as AGN 
feedback, known to have occurred in NGC 4649 dPaggi et alj2014h . 
can flatten the global metallicity p rofile because they quen ch star 
formation in the innermost regions dHirschmann et alj2014h . 

Given these continuous transformations of the metallicity pro¬ 
files, a direct comparison between the observed colour profiles of 
stars and GCs with galaxy simulations is more uncertain. Dry ma¬ 
jor mergers (our favoured formation model so far) can leave the 
metallicity of the remnant unchang ed if both progenitors h ave suf¬ 
ficiently steep metallicity gradients dPi Matteo et al.l2009bh . In this 
context, the observed negative gradients in NGC 4649 could be 
interpreted as the pre-merger gradients of the progenitor galaxies, 
which remained unchanged after the merger. This assumes that the 
merger triggered no star formation and that the initial metallicity 
slopes of the progenitors were equally steep. 

Our results show evidence of flat colour profiles out¬ 


side 20 kpc, consiste n t with other s t udies ([Rhodej&Zend 


Cantiello et al.l ^OTI: ^ssino et alj l2006l : Iporbes et ^ 


2001 


2011 


Faifer et ^ 201 iL Liu et al.l 2011 1. This feature is interpreted as 


the transition from the galaxy central part, created by dissipative 
processes, to an extended component made by secular accretion 
of (mostly metal-poor) satellite galaxies. In this context, GC bi- 
mod ality is a na t ural manifestation of dual galaxy formation path¬ 
way dC6tdll99^ : lBrodie et al.ll2014h : red GCs mostly form in-situ, 
along with bulk of stars, and are therefore more centrally concen¬ 
trated. Blue GCs are mostly accreted and are therefore more ex¬ 
tended ( lToninill2013h . 

The GC system of NGC 4649 is consist ent with being old 
(> 10 Gvr: lchies-Santos et alj201lh . However. IPierce et al]i2006ll 
found about four GC with very young ages in NGC 4649. These ob¬ 
jects are listed in Table |4] From this Table, we can see that young 
GCs have all red colours, but the sizes are ordinary for normal GCs. 
Given the lack of recent star formation in NGC 4649, the existence 
of these young GCs could be expla ined if they f ormed along with 
their host satellite galaxy at z ^ 1 ( lToninill2013h . which was then 
accreted onto NGC 4649. 

Lastly, we note that NGC 4649 was classified as a 
galaxy with “sub t le SO attributes” dSandage & TammannI Il98ll : 
ISandage & Bedkel Il994l). Processes such as secular evolution 
dLarson et al. 1980i r ^osellj_&Cjavazzil 200^; 
or gas stripping (e.g., iMoore et al.! Il996l: 


tipping 

11200711 


Kormendvl 

2013) 

Lake et alj 

1991 


IChung et^ 

a gas-poor SO i Kormendv & Bend^l201^ . 


can both transform a g as-rich spiral galaxy into 


ID - ID(P-|-06) 

Age 

[Gyr] 

V 

[km s“^] 

(9 - 2) 

[mag] 

rh 

[pc] 

R 

[kpc] 

J187- 1443 

2.1 ± 0.6 

703 ± 46 

1.53 

2.52 

3 

A8 - 502 

3.0 ±0.8 

705 ± 41 

1.50 

1.87 

7 

B36- 175 

2.4 ±0.4 

836 ± 11 

1.392 

1.95 

16 

B85 - 89 

2.4 ±0.8 

1148 ± 16 

1.29 

2.10 

18 


Tabl e 4. Properties of y oung GCs in NGC 4649. ID(P-|-06) and age are 
frornJPierceetalJ_ 2006h . The ID, colour (g — z) and size are from 
dStrader eTallTOli ). The velocity v is from this paper. Errors on (g — z) and 
Th are less than 10 per cent. R is the projected distance from NGC 4649. 


However, galaxies thought t o be faded spirals have the pho¬ 
tometric stellar disk still in place l| Cortesi et al.||2011|: [Forbes et al.l 


l2012l : lAmold et alj201ll : lcortesi et alj2013l : Kormendvll2013h . Nb 


strong evidence for a stellar disk has been found in NGC 4649 
dVika et alj 20 ITijKormendvl 20131) . Recently. I Graham et al.l d2015l) 
pointed out that NGC 4649 is particularly difficult to decom¬ 
pose into bulge and disk. It was recently shown that gas-rich ma- 
jor mergers can also p roduce SO remnan ts with embedded disks 
dOuereieta et al.ll2014h . iNaab et al.l d2014h predicts that these rem¬ 
nants have disky kinematics with rising or peaked Vrot jo profiles, 
along with hz and /14 strongly correlated with H/ct, which is not 
consistent with our findings. 


12.2 Final remarks and the interaction with NGC 4647 


We have shown that a formation pathway via a major merger be¬ 
tween two gas-poor galaxies, combined with satellite accretion, can 
simultaneously explain both the inner stellar observations and the 
outer GC observations. After a major-merger, the angular momen¬ 
tum is transported to the outer regions of the remnant, explaining 
the high rotation velocity observed for stars, GCs and PNe. 

If no gas was involved in the merger, the GC system of the 
primordial galaxies got mixed together, preserving their original 
metallicity gradients. Additional GCs (mostly metal-poor) are then 
accreted into the halo of the galaxy, explaining the extended spatial 
distribution of blue GCs. This scenario is supported by the presence 
of a depleted stellar core in the galaxy’s innermost regions. 

The possible interaction between NGC 4649 and NGC 4647, 
which is only 13 kpc away in projection, has been widely debated 
in the literature. A gri color image of this galaxy group is shown in 
Figure [TS] Radio observations suggest that the morphology of CO 
and HI maps are mildly disturbed, and so are the c orresponding ve¬ 
locity fields ( iRubin et all 19991 : 1 Young et al.l2006l) . However, radio 
observations alone do not return conclusive results on whether or 
not this gal axy pair is actually interacting . A similar conclusion was 
reached by Ide Griis & RobertsonI ( l2006l) from an analysis of ACS 
imaging of this galaxy. Moreover, the distribution of optical light 
around NGC 4647 and NGC 4649 shows no sign of disturbances 
(Figure [T^ . On tbe other hand, new observations found evidence 
for enhance d star formation i n tbe region betwe en NGC 4647 and 
NGC 4649 jLanz et aDl2013l:lMineo et al.ll2014t) . _ 

iD’Abrusco et al.l 2014h and Mineo et ^ ( 1201 4b have sug¬ 
gested that the substructures in the 2D distribution of GCs from 
ACS imaging may reflect a recent interaction in NGC 4649, but the 
role of NGC 4647 is still unclear. No significant sign of ongoing 
interaction between NGC 46 49 and NGC 4647 w as observed in the 
2D velocity field of the GCs dCoccato et al.l2013h based on the GC 
catalogue of iLee etal.l ( l2008h . 

Our kinematic results show some interesting features, such as 
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Figure 18. Subaru/Suprime-Cam gri colour image of NGC 4649 (left) and NGC 4647 (right). The projected separation between the two galaxies is 13 kpc. 
There are no clear signs of interaction between the two galaxies. North is up. East is left. 


wiggly velocity profiles, counter rotation, and a multi-spin PN sys¬ 
tem. We find that GCs with (g — z) ~ 0.9 ±0.1 mag, counter¬ 
rotate with respect to main rotation direction of the galaxy, whic h 
may repre sent the impr i nt of a recent interaction dSeth et al.ll2014h . 
Moreover, IPierce et"^ bOOfil) found a handful of GCs only 2 Gyr 
old in NGC 4649, which may be linked to a recent accretion 
event. Similar kinematic anomalies, such as strong rotation am¬ 
plitude and multi-spin features, have been found in GC systems 
of int eracting early-type galaxies, such as N GC 4365 l Blomet^ 
2012h. and merger remna nts like NGC 5128 jWoodlev et al.ll201^ 


Woodley & HarrisI 1201 ih . Computer simulations may determine 


whether the kinematic features observed in NGC 4649 encode 
signatures of an ongoing interaction between NGC 4649 and 
NGC 4647. Based on these contrasting results, it appears that, if 
NGC 4647 and NGC 4647 are interacting, we are witnessing the 
very early stages of this interaction. 


13 SUMMARY 

In this paper, we have studied the globular cluster (GC) sys¬ 
tem of the massive elliptical galaxy NGC 4649 (M60). We 
combined wide-field spectroscopy with wide-field imaging from 
space-based and ground-based telescopes. We build up a pho- 
tometric catalogue combining space-based HST/ACS data from 
I Strader et alj ( 1201 2h . ground-based archival CFHT/MegaCam, and 
new Subaru/Suprime-Cam data. GCs are selected based on colour 
and magnitude constrains. We detect all the photometric features 
usually seen in GC systems of large ellipticals: the GC system of 


NGC 4649 is made up of two subpopulations: the red (metal rich) 
GCs are centrally concentrated. Their spatial distribution is consis¬ 
tent with that of the starlight and planetary nebulae in NGC 4649. 
The blue (metal poor) GCs are more extended. The colour gradi¬ 
ents of blue and red GCs are flat outside 20 kpc, interpreted as 
evidence of accretion of metal-poor satellites into the stellar halo 
of NGC 4649. 

The spectroscopic follow-up is performed with three 
multi-object spectrographs: Keck/DEIMOS, Gemini/GMOS and 
MMT/Hectospec. These cover different effective areas around 
NCC 4649, but they are complementary with each other. We con¬ 
firmed 431 GCs associated with NGC 4649, plus some GCs from 
surrounding galaxies, and dwarf galaxies. The contamination from 
stars and background galaxies is consistent with zero. We analyse 
the kinematics of blue (60 per cent of the total sample) and red GCs 
(40 per cent), along with literature planetary nebulae data. Results 
are: 

• Significant rotation at all radii is detected for stars, PNe and blue 
and red GCs. Rotation occurs preferentially along the major axis of 
the galaxy. 

• The velocity dispersion of the three tracers look very different. 
Blue GCs and PNe have higher velocity dispersion in the outer re¬ 
gions, whereas the red GCs and the innermost PNe are consistent 
with the stellar kinematics and have lower dispersions. 

• We report peculiar kinematic features for both GCs and PNe. The 
rotation axis of PNe twists with radius, whereas a group of blue 
GCs at intermediate radii counter-rotate with respect to the bulk 
of blue GCs. The velocity dispersion profiles are bumpy, a feature 
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which may indicate that this galaxy underwent major mergers in 
the past. 

We discuss possible formation scenarios for this galaxy, con¬ 
trasting stellar and GC observations with recent galaxy formation 
simulations. We find that a formation via dry (no gas) major-merger 
between two galaxies can, regardless of their initial angular mo¬ 
mentum, consistently reproduce stellar observations in the inner¬ 
most regions and GC observations in the outermost regions of 
NGC 4649. We consider the possibility that NGC 4649 was spi¬ 
ral galaxy which evolved in an elliptical galaxy because of gas 
stripping, but the apparent non-detection of the primordial disk in 
NGC 4649 makes this scenario unlikely. We find no strong evi¬ 
dence to support an upcoming interacting between NGC 4649 and 
the spiral NGC 4647. 
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